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A detailed study of the slow manifold of a model exhibiting mixed-mode oscillations is presented. 
A scenario for the emergence of mixed-mode states which does not involve phase locking on a 
2-torus is constructed. We show that mixed-modes correspond to the periodic orbits embedded 
in the horseshoe-type strange set and demonstrate how transformations of this set determine the 
transitions of mixed-mode states into each other. 


I. INTRODUCTION 

Oscillatory behavior exhibited by nonlinear chemical 
systems often takes the form of periodic mixed-mode 
oscillations consisting of L large amplitude oscillations 
followed by S small amplitude oscillations. The distinc¬ 
tive difference between large and small oscillations allows 
one to classify the mixed-mode oscillatory states by in¬ 
teger pairs ( L,S ) and assign them the symbol L s . Or¬ 
dered progressions of such periodic patterns separated 
by chaotic oscillations are observed as control param¬ 
eters are varied. The progressions of different mixed¬ 
mode forms are known to be Farey sequences described 
by Farey arithmetic [lj] . 

Since their discovery in the Belousov-Zhabotinsky re¬ 
action [2j, mixed-mode oscillations have been found in 
a broad range of chemical and biological [|J systems. 
Their ubiquity suggests the possibility of the existence of 
a common dynamical systems theory origin. A number 
of authors |lJU identified experimentally observed mixed¬ 
mode oscillations with periodic motions arising as a result 
of phase-locking on a 2-torus. This hypothesis was based 
mainly on the fact that the periodic states on an invari¬ 
ant torus form Farey sequences as a control parameter 
runs through the quasiperiodic route from period-1 oscil¬ 
lations to chaos. Numerical evidence for the existence of 
quasiperiodicity and torus oscillations was obtained from 
studies of model systems which exhibit mixed-mode os¬ 
cillations Hi). 

However, there are results which show that mixed¬ 
mode oscillations cannot be attributed to the phe¬ 
nomenon of quasiperiodicity in all systems. In a study of 
electrodissolution of copper Albahadily et al. JtJ demon¬ 
strated that the mixed-mode states emerge in parame¬ 
ter regions different from those where a torus is stable. 
The formation of mixed-mode oscillations was observed 
at a point well beyond the range of parameters between 
the secondary Hopf bifurcation and the torus break-up 
bifurcation, where the mixed-mode states should lie ac¬ 
cording to the hypothesis of quasiperiodicity. In a series 
of electrochemical studies Koper et al. ||] found that the 


mixed-mode oscillations in Farey sequences observed in 
their study are separated by chaotic states which resem¬ 
ble random mixtures of the adjacent periodic patterns 
rather than by quasiperiodic oscillations densely cover¬ 
ing the surface of a 2-torus. In a numerical study of a 
model system Koper |^] found a parameter domain where 
both quasiperiodic oscillations lying on a 2-torus and a 
large-amplitude mixed-mode oscillatory state which does 
not belong to the torus are stable simultaneously. In or¬ 
der to interpret the results of the experimental study in 
Ref. [Q Ringland et al. |l(J constructed a family of two- 
extrenrum Z-maps capable of producing Farey sequences 
in the limiting case when the central segment of the Z- 
map is vertical. They have shown that one may obtain 
Farey sequences which are not related to phase locking 
on a torus. 


In this paper we present results of investigations of 
a model proposed earlier for the qualitative description 
of the mixed-mode oscillations observed in the Belousov- 
Zhabotinsky reaction |ll|]. From the analysis of the trans¬ 
formation of the slow manifold of the model as system pa¬ 
rameters are changed we demonstrate the gradual emer¬ 
gence of a horseshoe-type strange set in which all the 
mixed-mode states are embedded as periodic orbits. On 
the basis of this analysis we construct a three-dimensional 
phase space picture describing the emergence and bifur¬ 
cations of mixed-mode states. This scenario does not 
involve a torus and appears to be an alternative to the 
well-established quasiperiodic scenario. 


In Sec. [0] we introduce the model and present a brief 
discussion of its attractors and bifurcation structure. 
Section |l| is devoted to the detailed analysis of the slow 
manifold and a description of its important properties. 
We introduce a technique that allows us to construct the 
slow manifold and demonstrate how it turns into a horse¬ 
shoe by stretching and folding . The repulsive flow that 
defines a particular shape of mixed-mode oscillations is 
introduced and discussed in some detail. In Sec. IV we 
show that when the slow manifold is folded into a horse¬ 
shoe interacting with the repulsive flow, the mixed-mode 
oscillations arise naturally at the point where the large- 
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amplitude, period -1 oscillation loses its stability through 
a saddle-node bifurcation. The relation of the phase flow 
on the horseshoe to the Z-map introduced by Ringland 
et al. [0 is discussed. In Sec. |v| we conclude that the 
existence of the repulsive flow allows one to consider this 
horseshoe as a representative of a distinct subclass of 
strange sets whose periodic windows are organised ac¬ 
cording to Farey sequences. 


II. DESCRIPTION OF THE MODEL 


We consider a three-variable extension El of the 
Boissonade-De Kepper model |h|, 

x = 7 (y — x 3 + 3yx) , 
y = -2yx-y-z + (3, (1) 

z = k(x — z ) . 

System (§ was originally proposed in | jTl| in order to 
obtain a qualitative description of transient mixed-mode 
oscillations in the Belousov-Zhabotinsky reaction. De¬ 
pending on the values of y and (3, (Jli) possesses either 
one or three steady states. The change in the number of 
stationary points occurs on the curve 



where two steady states coalesce or emerge in a saddle- 
node bifurcation. At y = 1, (3 = 0, where two branches 
of the curve (|j) meet in a cusp point, the pitch-fork bi¬ 
furcation takes place. Note that the parameters 7 and 
ft only rescale the corresponding components of the vec¬ 
tor field and do not influence the positions of the null- 
clines and steady states in phase space. In the following 
discussion the geometric parameters /i and (3 are fixed 
(y = 2.0, (3 = —0.4) while 7 and ft are varied. With 
this choice of y and [3 system (jlj) has only one steady 
state at Xf = Zf = —1.159705, yf = 5.398524, which is 
a stable focus in a domain of the parameter plane ( 7 , ft) 
lying to the left of the Hopf bifurcation line (see Fig. [I]), 
and a saddle-focus otherwise. As the real characteris¬ 
tic exponent of the steady state is negative everywhere 
on ( 7 , ft), no second order bifurcation is possible at the 
chosen values of y and f3. 

Figure [I] shows the partition of the ( 7 , ft) plane into 
domains of different types of dynamical behavior accord¬ 
ing to the results of numerical integration. In domain 
1 a stable focus is the only attractor of the system 0 )- 
Large amplitude period-1 ( 1 °) oscillations with a strong 
relaxational character at large 7 are found in domain 2 . 



FIG. 1. Bifurcation phase diagram: 1 - steady state (fo¬ 

cus), 2 — period-1, 3 — mixed-mode oscillations. Separat¬ 
ing lines correspond to the Hopf bifurcation (dots), transition 
1 ° <-> oo 1 (dashes), and crash onto the focus (solid). Both 
steady state and oscillations are stable in subdomains 2 ' and 
3'. 

Finally, mixed-mode oscillations are observed in the 
wedge-shaped domain 3. Transitions from period-1 to 
mixed-modes take place on an almost flat curve (shown 
by the dashed line in Fig. [I]) where the trajectory of sys¬ 
tem (|l|) performs a small-amplitude excursion only once 
on the entire time interval t G (— 00 ,+ 00 ). This transi¬ 
tion can be represented symbolically as 1 ° <-» oo 1 . Con¬ 
tinuing to traverse domain 3 by decreasing 7 at constant 
ft one observes pruned Farey sequences of periodic mixed¬ 
mode oscillations oo 1 —> n 1 —> l 1 —> l n —> I 00 inter¬ 
spersed by narrow, yet detectable, spans of chaos. This 
periodic-chaotic progression ends suddenly in a crash 
onto the focus at the point where the line separating 
zones 1 and 3 is crossed. At large ft (ft > 0.54), where 
the crash bifurcation curve appears to be tangental to 
that for the 1 ° <-> oo 1 transition, both bifurcations oc¬ 
cur within an extremely narrow parameter span and one 
observes hard termination/emergence of large-amplitude 
periodic oscillations. Another limiting case is observed at 
low ft (ft < 0 . 20 ) where the crash curve merges with the 
Hopf bifurcation line. Here the crash occurs from pattern 
l n , n>l which forms an almost homoclinic connection 
to the focus undergoing the Hopf bifurcation. 

As ft varies the appearance of the phase portrait of 
the mixed-mode states changes considerably. At large 
values of k both small- and large-amplitude excursions 
lie almost on the same plane. This type of mixed-mode 
oscillation, referred to type -1 according to classification 
given in ||], is shown in Fig. ||(a,b). With a decrease in ft 
the reinjection part of trajectory falls closer to the focus 
and the angle between the planes of small- and large- 
amplitude oscillations grows. 
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FIG. 2. Projections of the phase portraits on ( x , z ) 
plane and corresponding x(t) time series for two periodic 
mixed-mode states: (a,b) — 2 1 for k = 0 . 5,7 = 0.35516; 
(c,d) - l 12 for k = 0.2 ,7 = 0.471. 

At k < 0.2 the phase portrait of the oscillation acquires 
all of the features characteristic of type -2 mixed-mode 
states. Here the large loops are mainly orthogonal to 
the small ones and the reinjection occurs along the one- 
dimensional stable manifold of the focus (see Fig. @(c,d)). 

III. CONSTRUCTION OF THE PHASE FLOW 

The parameter region in which complex periodic and 
chaotic oscillations are found can be partitioned into 
subdomains corresponding to particular periodic mixed¬ 
mode states L s and their chaotic mixtures. We shall 
study the mechanisms responsible for the bifurcations of 
one mixed-mode state into another by analysing transfor¬ 
mations of the phase flow on the system’s slow manifold. 

The concept of a slow manifold relies on the assump¬ 
tion that the relaxation of an arbitrary initial condition 
to the attractor has two characteristic time scales: the 
trajectory quickly reaches certain subset of the phase 
space and then slowly approaches the attractor within 
this subset. In the case of a strongly dissipative system 
with three dynamical variables, the slow set is usually 
represented by a 2D manifold which is often folded in 
the 3D phase space. To avoid complications arising from 
nontrivial embedding of the slow manifold, we choose for 
its construction the region of large k where the slow man¬ 
ifold appears to have a relatively simple organization. 

A. Slow manifold 

Figure ||(a) presents the projection of the phase flow 
on the (x, y) plane for n = 0.55 and 7 = 0.39. For these 
parameter values deep within domain 2 ' the system ( 0 ) is 
bistable. Apart from the stable focus, two stationary pe¬ 
riodic solutions exist: a stable, large-amplitude, period -1 
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FIG. 3. Phase flow on the slow manifold, (a) Projection 
of the flow on (x , y) plane, (b) Mutual positions of the slow 
manifold and the surface-of-section P in phase space, (c) First 
return map for the flow initiated on Aq. 
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limit cycle (solid line) and a saddle-type, small- 
amplitude cycle (dashed line). The planar, concentric 
organization of these cycles around the focus suggests 
that the slow manifold containing them is flat. Trajecto¬ 
ries that leave the saddle cycle at t — > — oo and approach 
stable cycle or focus at t — > +oo lie in the slow man¬ 
ifold. These trajectories are easy to calculate and one 
can use them for the construction of a Poincare section 
of the slow manifold. We choose P = {x = Xf, y < yf} 
shown in Fig ||(a,b) as a surface of section. Intersec¬ 
tions of the transient trajectories with P yield a set of 
points which, after interpolation by a smooth function 
z = z(y) |l4j], becomes a continuous zero-order approxi¬ 
mation Aq = {x = Xf 7 y < yf, z = z(y)} to the Poincare 
section of the slow manifold. 

Figure ||(b) shows a 3D phase portrait of the same 
trajectories as in Fig. ||(a), as well as P and Aq. The in¬ 
tersections of the flow initiated on Aq with P can be used 
both to improve the accuracy of the approximation and 
to construct the first return map (see Fig. ||(c)) which 
quantitatively describes the phase flow on the slow man¬ 
ifold. 

A complete description of the phase flow on the slow 
manifold can be obtained from the comparison of the 
phase portrait and the first return map. Consider 
the example presented in Fig. ||(a,c). The first re¬ 
turn map y n + 1 = f(y n ) intersects the bisectrix in three 
points which correspond to the stable limit cycle (y c = 
—2.5437), the saddle cycle {y s = 4.2998), and the stable 
focus (yf = 5.3985). The nearly rectangular shape of 
the map implies that the slow manifold comprises areas 
of strong expansion (the neighborhood of yo = 3.9370, 
where the map is vertical) and contraction (horizontal 
part of the map). In the course of its evolution a trajec¬ 
tory initiated in the vicinity of the saddle cycle executes 
a number of small-amplitude loops around the saddle cy¬ 
cle and then in one turn reaches the neighborhood of the 
stable limit cycle. Whether the trajectory approaches 
the limit cycle like trajectory 3 performing small 



FIG. 4. Scaling of f(y) at the point yo where tangent is 
vertical. 


transient loop or like trajectory 4 which makes a large 
loop depends on its position relative to a certain bound¬ 
ary trajectory 0 which can be identified with the point 
yo on the first return map where f(y) has a vertical tan¬ 
gent. Analysis of f(y) in a small neighborhood of yo 
shows that the trajectory 0 is superunstable; i.e., the 
neighboring trajectories deviate from it faster than geo¬ 
metrically. Figure ^ shows log 10 |A/| versus log 10 \Ay\, 
where A/ = f(y) - f(y 0 ) and Ay = y - y 0 , for both 
left y < yo (lower set of points) and right y > yo (upper 
set) branches of the map shown in Fig. 0(c). One can see 
that for Ay sufficiently small, A/ (Ay) is characterized 
by the power law A/ oc (A y) a where a « 0.022 <C 1 . 
This implies an extreme divergence of trajectories which 
is also reflected in the phase portrait. The initial con¬ 
ditions (x f ,yi,z(y 1 )), (xf,y 2 ,z(y 2 )) on P for trajecto¬ 
ries 1 and 2 (solid lines in Fig. |3|(a,b)) were such that 
= 2 /o - S/2, 2/2 = 2 /o + S/2 where <5 = 10 -12 . As Fig. § 
shows, the separation between these two trajectories on 
their return to the surface of section grows by « 12 or¬ 
ders of magnitude. In fact, trajectories 3 and 4 constitute 
another pair with S = 10~ 10 . 

The existence of the superunstable trajectory results 
in the observed dichotomy of trajectories belonging to 
the slow manifold. A 2D manifold intersecting the slow 
manifold along trajectory 0 plays similar separating role 
in the 3D phase space. Although the global structure 
of this repulsive flow in the phase space is very compli¬ 
cated, the flow can be easily determined locally by the 
property that it separates trajectories making large and 
small loops while approaching the attractor. 


B. Birth of horseshoe 

Figure ||(a) shows an expanded view of the flat part 
of the first return map f(y) in Fig. ||(c). The slope 
of the map is everywhere positive and, thus, the map 
is invertible. As 7 decreases at constant k = 0.55, a 
smooth transformation of the slow manifold takes place, 
resulting in the loss of invertability of the map. A 
fragment of map constructed for 7 = 0.332241 (see 
Fig. |](b)) demonstrates that f(y) is not monotonic and 
a smooth maximum arises near the stable limit cycle. At 
7 = 0.33224025 the map is tangent to the bisectrix and 
a new pair of stable and unstable limit cycles emerges on 
the slow manifold (cf. Fig. ^j(c)). As the map’s maximum 
grows it sharpens and the newborn limit cycle looses its 
stability through a cascade of period-doubling bifurca¬ 
tions. At 7 = 0.3322392 (see Fig. ||(d)) a chaotic at¬ 
tractor coexists with the period -1 limit cycle and focus. 
Although the entire cascade to chaos occupies very nar¬ 
row span of 7 , two first period-doubling bifurcations were 
resolved. The tristability of the system breaks down at a 
certain value of 7 and every trajectory initiated near the 
maximum of the map is attracted to the period -1 limit 
cycle after a chaotic transient of varying length. 
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FIG. 5. Transformation of the first return map at k = 0.55 
and decreasing 7 : (a) 7 = 0.39, (b) 7 = 0.332241, (c) 
7 = 0.3322399, (d) 7 = 0.3322392, (e) 7 = 0.332239127. 


Figure ^](e) shows a part of the first return map for 
7 C = 0.332239127 where oscillations suddenly crash on 
the focus. Two important characteristic features of the 
flow can be inferred from the map. First, the maximum 
of the map acquires a cuspoid shape signalling of the for¬ 
mation of a highly compressed fold on the slow manifold. 
Second, from the tangency of the map and bisectrix, one 
can see that the stable period -1 limit cycle and the large- 
amplitude saddle cycle are near the point of coalescence. 
The annihilation of this pair in a saddle-node bifurcation 
gives rise to an intermittent chaotic attractor which has 
the appearance of a thick limit cycle. A decrease in the 
control parameter 7 leads to the crash of chaotic oscilla¬ 
tions onto the focus. Summarising the information pre¬ 
sented in Fig. ||(a-e), one can conclude that as 7 changes 
through the interval [0.39, 7 J, the slow manifold trans¬ 
forms from a planar surface into a folded fractal structure 
with infinitely many leaves (Smale’s horseshoe). 

Making use of several auxiliary Poincare sections it 
is possible to construct a 3D picture of the slow mani¬ 
fold. Figure |] schematically represents the manifold and 
three of its Poincare sections at parameters close to the 
crash bifurcation. Part of the manifold situated between 
sections C and A is removed to facilitate visualization 
of its geometry. In section A one sees a triply-folded 
curve which corresponds to the first structural level of 
the manifold. Note the difference between the two folds 
of this curve: while the outer fold has a smooth parabolic 
shape, the inner fold is a sharp cusp. Increasing the res¬ 
olution considerably (by « 5 orders of magnitude) one 
would be able to find that the curve is, in fact, a densely 
compressed fold constituted by seven leaves of the sec¬ 
ond structural level, and so on. To describe one iteration 
of the manifold stretching and folding we do not need 
to consider these higher structural levels. Following the 
evolution of points in section A under the flow (clockwise 
in Fig. |g]) one encounters a pleat where the manifold as a 
whole begins to fold into an ^-shaped form (section B) 
with two approximately equal parabolic folds. Prominent 
stretching of the flow in the plane of the manifold, 



FIG. 6 . Geometry of the folded slow manifold. The supe¬ 
runstable trajectory is shown by dot-and-dash line. 
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accompanied by strong contraction in the orthogonal 
direction, transforms the manifold into the thin roll com¬ 
posed of seven leaves seen in section C. As a result of 
the extreme inhomogeneity of compression the inner fold 
acquires sharp cuspoid shape while the outer one remains 
smooth and parabolic. As the flow returns to section A 
more contraction along the vertical direction occurs and 
the second order structure shown in section C is not dis¬ 
cernible. 


IV. MIXED-MODE OSCILLATIONS 


In Sec. [m] we described two important characteristic 
features of the phase flow of system (Lj). Firstly, at the 
parameter values in the neighborhood of mixed-mode os¬ 
cillations, the slow manifold is folded into a horseshoe- 
type set with a fractal, multi-leaved structure. Secondly, 
the repulsive flow divides the transient trajectories into 
two distinct classes. Intersection of the repulsive flow 
and the slow manifold yields a superunstable trajectory 
with peculiar dynamical properties. In the present sec¬ 
tion we show how mixed-mode states which correspond 
to periodic orbits embedded in the horseshoe emerge and 
bifurcate into each other. 


A. Emergence of mixed-mode oscillations 


Consider a path in the ( 7 , k ) parameter plane which 
starts in domain 2 or 2 ' and ends in the region of mixed¬ 
mode oscillations 3 (see Fig. 0). As the domain of large- 
amplitude period -1 oscillations is traversed towards the 
curve corresponding to the 1 ° <-> oo 1 transition, the mod¬ 
ifications of the slow manifold described in Sec. Ill B take 
place. Before the first mixed-mode is formed an impor¬ 
tant transformation of the manifold occurs at the point 
where a cusp-shaped fold developes. Inspection of Fig. |0 
shows that the manifold cannot support periodic oscil¬ 
lations with small-amplitude loops since its folded part 
lies on the “large-amplitude” side of the repulsive flow 
(only the superunstable trajectory is shown in Fig. 0). 
In order for orbits with both small- and large-amplitude 
excursions to exist, the repulsive flow has to intersect the 
folded part of the horseshoe. 

Figure 0 shows segments of the manifolds cut out by 
two orthogonal planes for three consequtive values of 7 
and fixed k. The case 7 > 7 c (Fig. 0(a)) corresponds to 
the generic transversal position of the manifolds in which 
the growing fold of the horseshoe lies on one side of the 
repulsive flow. This case is described by a first return 
map of the type shown in Fig. 0(d). At 7 = in addi¬ 
tion to the transversal intersection, a tangency between 
the two manifolds occurs. The tangency manifests itself 
by the formation of a sharp cusp on the first return map 
(see Fig. 0(e)). For 7 < 7 C the fold of the horseshoe pen¬ 
etrates the surface corresponding to the repulsive flow 




FIG. 7. Propagation of the slow manifold through the sur¬ 
face of the repulsive flow (hatched): (a) 7 > -y c — transversal 
intersection; (b) 7 = 7 C — formation of a tangency; (c) 7 < 7 C 
— fold penetrates the surface. 


and drastically increases in size due to the strong diver¬ 
gence from the repulsive flow. One also sees in Fig. 0(c) 
that as a result of this transformation two new supe¬ 
runstable orbits (thick solid lines) are born on the slow 
manifold. Taking into the consideration fractal nature 
of the horseshoe one could imagine an infinite stack of 
superunstable orbits born in pairs when nested folds of 
the horseshoe penetrate the surface of the repulsive flow 
one by one. 

The birth of new superunstable orbits is reflected in the 
first return map as well. Figure |0 shows the first return 
map construct ed wi th P as a surface of section as de¬ 
scribed in Sec. Ill A for parameter values k = 0.50, 70 = 
.355403274 < 7 C at which the transition 1 ° <-> oo 1 occurs. 
As 7 C is passed, the cusp on the map breaks up and the 
maximum acquires normal parabolic shape. At the same 
time the height of the maximum increases tremendously 
and two segments with infinite slope appear on its sides. 

The analysis of the map shows that the transition 
1 ° <-► oo 1 occurs simultaneously with annihilation of the 
large-amplitude stable and saddle-type limit cycles. At 
the point of coalescence they produce structurally un¬ 
stable limit cycle which is seen on the first return map 
as a point of tangency with bisectrix. This cycle com¬ 
bines properties of both stable and unstable limit cycles. 
Consider a trajectory which departs from such cycle at 
t —> — 00 . In course of evolution it inevitably executes 
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FIG. 8 . First return map at the point of the 1 ° <-> oo 1 
transition (k = 0.5, 7 = 0.355403274). 

one small-amplitude loop around the focus. Upon the 
completion of this loop the trajectory is reinjected into 
the stable part of the neighborhood of the cycle which is 
eventually approached at t —> + 00 . Thus, first mixed¬ 
mode state oo 1 is characterized by the orbit homoclinic 
to the structurally unstable cycle and, therefore, exists 
in a parameter domain of zero measure. For 7 < 70 a 
finite number of large-amplitude loops is executed before 
a trajectory visits the small loop and successions of n 1 
mixed-mode states with rapidly decreasing n are oserved. 

B. Embedding of the mixed-mode orbits 

It is possible to establish how a particular mixed-mode 
orbit is incorporated in the slow manifold. Figure shows 
the (x,y) projection of a 46 1 orbit and its embedding 
in the horseshoe. As in Fig. |(| the Poincare sections 
of the manifold are shown at different levels of resolu¬ 
tion of the fine structure. While section A presents the 
first two structural levels, in section B only first level is 
shown in full. Two lines of phase points parallel to the 
surface of the repulsive flow and separated from it by 
S = ±0.5 x 10 ~ 12 are selected in section A (long dashes). 
The evolution of these points shown in section B demon¬ 
strates the extreme divergence of trajectories caused by 
the presence of the repulsive flow. 

Orbit 46 1 consists of a dense, flat band of large- 
amplitude loops which consecutively shrink in size and a 
small loop where the trajectory segment is reinjected into 
the large-amplitude band. Projecting points in which the 
orbit intersects planes A and B (shown by thick dots) on 
the corresponding Poincare sections of the horseshoe one 
can obtain insight into the mechanism of this reinjec¬ 
tion. At a certain time the trajectory appears to the left 
of the repulsive flow surface (point 1 in A). This causes 
it to execute a small loop (point 2 on the cusp-shaped 
tongue in B). Next return of the trajectory into section 
A occurs in point 3, separated far from 1. Following the 


subsequent numbers one sees how the orbit returns to the 
large-amplitude band, gradually ascending leaves of the 
slow manifold. The same mechanism holds for all n 1 




FIG. 9. Projection of 46 1 orbit onto (x , y) plane and its 
embedding in the slow manifold (see discussion in the text). 
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states, although escape on the small loop occurs from 
different leaves. 

The analysis of the reinjection scheme shows that the 
horseshoe (solid lines in Fig. |^) cannot support mixed¬ 
mode oscillations with more than one small-amplitude 
loop. Indeed, regardless of how far to the left one shifts 
the position of point 1 , first return of a trajectory initi¬ 
ated at 1 into the surface of section (point 3) lies to the 
right of the repulsive flow. This situation changes as the 
cusp-shaped tongue, to which point 3 belongs, penetrates 
the surface of the repulsive flow t hroug h the formation 
of a tangency as described in Sec. IV A (Corresponding 
changes on both sections are shown by medium-dashed 
lines.) As soon as a similar tongue of a higher struc¬ 
tural level passes through the surface of the repulsive 
flow, an opportunity for the existence of the third small- 
amplitude loop arises, and so on. Thus, changing po¬ 
sition of the horseshoe relative to the repulsive flow one 
can obtain mixed-mode states L s with any desired L and 
S. 


C. Relation to Z-map 


At small values of k, where the bulk of the mixed¬ 
mode domain lies, the surface-of-section P introduced in 


Sec. Ill A becomes unsuitable for the construction of the 


first return map due to the strongly bent shape of the 
corresponding Poincare section (see section B in Fig. |^). 
Instead, one can use sections y = const similar to sec¬ 
tion A in Fig. ||. Figure [To] presents two examples of the 
first return map x n +i = fz{x n ) constructed in the ap¬ 
proximation based on the upper leaf of the horseshoe for 
parameter values corresponding to 46 1 (a) and l 1 (b) 
mixed-mode states. 

As one can see, maps of this type share a number of 
common features. The map consists of two branches with 
positive slope and an extremely steep, negatively-inclined 
segment which joins the branches. At the point corre¬ 
sponding to the superunstable orbit the map has an in¬ 
finite slope and can b e loc ally described by an exponent 
as discussed in Sec. Ill A . Any mixed-mode state L s 
may be represented on this map by a periodic trajectory 
with L iterates on the right branch and S iterates on the 
left branch. This property of the map makes it particu¬ 
larly suitable for the illustration of mixed-mode bifurca¬ 
tions. At parameter values corresponding to the transi¬ 
tion 1 ° <-» oo 1 the right branch is tangent to the bisectrix. 
As an intermittent channel opens, the trajectory funnels 
through it and then jumps onto the left branch from 
which it is reinjected onto the branch of large-amplitude 
loops and the cycle closes. As 7 decreases both branches 
fall relative to the bisectrix, reflecting the gradual growth 
of the number of small-amplitude loops accompanied by 
a reduction in the number of large-amplitude loops. 
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FIG. 10. First return maps fz(x ) constructed for k = 0.4 
with y = 0.4 as a surface-of-section: (a) 7 = 0.46204, (b) 
7 = 0.44. Orbits 46 1 and l 1 are shown with thin solid lines. 


In Ref. [|lOj Ringland et. al. extensively studied the 
properties of the two-extremum map = Z(x n ) 

where 

Z( x) = (c + ax){ 1 — tanh(sx)) + (d + bx)(l + tanh(sa;)). 

At small values of the parameter s the map has a smooth 
shape similar to that of the cubic map. As s increases, 
the slope of its middle, negatively-inclined segment grows 
steeper while the x distance between two map’s extrema 
vanishes. As s tends to infinity, the Z-map acquires 
zigzag shape with a vertical middle segment. Ringland 
et. al. have shown that in the limiting case s —> 00 the 
attractors of the Z-map form Farey sequences as other 
map parameters are varied. They related this property 
to the existence of a vertical segment since a decrease in 
s leads to a gradual transformation of Farey sequences 
into U-sequences. Although the dynamic origin of such 
maps was not discussed, a possible relation of the Z-map 
to the existence of mixed-mode oscillations was inferred. 

One can easily see that the Z-map in its s —> 00 limit 
is qualitatively similar to the map fz(x) constructed in 
the present section. In our map the middle segment is 
always vertical due to the presence of the superunstable 
orbit. Applying results of Ringland et. al. to system ( 1 ), 
we may assume that the repulsive flow defines not only 
the shape of mixed-mode orbits with their nonoverlap¬ 
ping bands of small- and large-amplitude loops, but also 
their property to form Farey sequences. 
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V. CONCLUSIONS 


VI. ACKNOWLEDGEMENTS 


We have analyzed in detail a model exhibiting mixed¬ 
mode oscillations. By constructing the model’s slow man¬ 
ifold we have shown that mixed-modes correspond to pe¬ 
riodic orbits embedded in a horseshoe-type strange set. 
This explains why chaotic oscillations are observed in 
transitions between adjacent periodic states. The organi¬ 
zation of the slow manifold into a horseshoe also accounts 
for the signatures of chaos in transient trajectories when 
the attractor is periodic. 

In our model the mixed-mode periodic orbits do not 
lie on a 2-torus. Following transformations of the slow 
manifold from its simple planar organization into a horse¬ 
shoe, we have also shown that a torus does not exist as an 
intermediate state between period -1 and mixed-mode os¬ 
cillations. Thus, our scenario for the emergence of mixed¬ 
mode oscillations is an alternative to quasiperiodicity. 

The main distinctive feature of mixed-mode states, the 
partition of their orbits into nonoverlapping bands of 
small- and large-amplitude loops, finds its explanation 
in the existence of the repulsive flow. Intersection of this 
flow with the slow manifold yields a complex system of 
superunstable orbits. The position of the slow manifold 
relative to the repulsive flow determines whether the for¬ 
mation of small or large loops of the attractor is favored. 
As parameters change along a path running through the 
mixed-mode domain, the position of the repulsive flow 
changes and one observes sequences of mixed-mode states 
n 1 —> l 1 —* l n (l n —> l 1 —>■ n 1 ) with monotonously in¬ 
creasing (decreasing) S/L ratio. As the results of Ring- 
land et.al. suggest, the fact that these sequences are 
described by Farey arithmetic is also a consequence of 
the presence of the superunstable orbits. Thus, the ex¬ 
istence of the repulsive flow places this horseshoe into a 
separate subclass of strange sets whose periodic orbits 
are given by L s mixed-mode states, and whose windows 
of periodicity are organized into Farey sequences. 

There exist indications that our scenario may account 
for the mixed-mode oscillations observed in a number 
of model and experimental studies. This appears to be 
so for all those cases in which the mixed-mode periodic 
states were found to be separated by chaotic rather then 
quasiperiodic oscillations. As Ringland et.al. jl(J pointed 
out, zigzag-shaped maps with a vertical middle segment 
indicative of a superunstable trajectory were found in 
many chemical and electrochemical systems. 


This work was supported in part by a grant from the 
Natural Sciences and Engineering Research Council of 
Canada. 
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